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Abstract. Two continuum subtracted Ha HST frames 
of M101 are used to determine the positions, angular 
sizes and absolute fluxes of 237 H II regions using a semi- 
automated technique. From these we have constructed the 
luminosity and diameter distribution functions. 

We repeat this process on the images after artificially 
reducing the linear resolution to that typically obtained 
with ground based imaging. We find substantial differ- 
ences in the luminosity function and diameter distribu- 
tion. The measured internal properties, such as central 
surface brightness and radial gradient are dominated by 
the PSF at linear resolutions less than roughly 40 pc 
FWHM. From the ground such resolutions are currently 
only obtainable for the nearest galaxies. 

Further support for the dominant role played by the 
seeing is provided by simple analytical models. We also 
study the clustering properties of H n regions and their 
effect on the luminosity function by construction of a 
Minimal Spanning Tree (MST). We find evidence for two 
regimes of clustering of the H II regions and diffuse emis- 
sion. These intrinsic clustering properties in combina- 
tion with the spatial resolution typically obtainable from 
ground based observations might be responsible for the 
break in the H n region luminosity function which is usu- 
ally found at logLna = 38.6erg/s, suggesting two different 
regimes of star formation in late type spiral galaxies. 

From the high resolution HST data we find a lumi- 
nosity function slope of a = —1.74 ± 0.08. We also ob- 
serve a flattening at luminosities logL < 36.7erg/s. For 
the diameter distribution we find a characteristic scale of 
i?o = 29.2 pc from an exponential fit. However, a scale free 
power law with index (3 = — 2.84± 0.16 provides a better 
fit to the data. 
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1. Introduction 

As tracers of recent massive star formation, Hll regions 
are indispensable. The radiation emitted by the hot gas 
carries the signature of the young hot stars to which 
it gave birth. Complementing the detailed spectroscopic 
work that has been done, are studies on the statistical 
properties of Hll regions. These global investigations are 
aimed at gaining an understanding of the connection be- 
tween the large scale dynamics and the star formation 
properties of galaxies (e.g. interactions and mergers, bars 
and spirals, resonances and rings). Most of this work has 
focused on the luminosity and diameter distribution of the 
regions imaged over full galactic disks in the line of Ha. 

A long standing difficulty with making such measure- 
ments has been in the identification of the regions. At low 
resolutions the blending together of regions can be sig- 
nificant. An initial estimate of this effect was made by 
Kennicutt et al. ( 1989] ), who found that blending effects 
became important at linear resolutions less than 300 pc. 
However, at the luminosity limits now obtainable with 
modern CCDs, a study of resolution effects has not been 
made. 

For the diameter distribution of Hn regions Van den 
Bergh (1981) found the functional form of an exponential 
law. Hodge (1987) confirmed this form for a sample of spi- 
ral galaxies and derived a dependence of characteristic size 
Dq on absolute magnitude of the galaxy. A wide range of 
of Do between 45 and 560 pc was observed for the different 
galaxies. Nevertheless, in other studies a power law form 
fit the diameter distribution better than the exponential 
function (Kennicutt & Hodge 



198C, Elmegrcen & Salzer 



M100 that is 1.5 to 2.2 times smaller than found in previ- 
ous determinations at lower resolutions. The reliability of 
available diameter measurements is highly questionable. 

A systematic study of the luminosity function (LF) of 
the H II re gions i n disk galaxies was carried out by Kenni- 
cutt et al. ( 1989 ), using a sample of 30 spiral and irregular 
galaxies. They found the LF could be represented by a 
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power law with index a — —2 ± 0.5. A flattening at lumi- 
nosities below logLna = 37erg/s was found and explained 
in terms of a transition from single ionizing stars to ion- 



izing stellar clusters (McKee & Williams 1997 and Oey & 
Clarke 1998). The LF slope was also found to be system- 



atically steeper for earlier Hubble types and the properties 
of regions in spiral arms differed to those between arms. 
Theoretical explanations for these variations were given 
by Oey & Clarke ( 1998] ), who showed that evolutionary 
effects and maximum number of ionizing stars per cluster 
can account for such differences. Some galaxies also exhib- 
ited a break in the LF at about logLna = 39erg/s with a 
steeper slope at the high luminosity end. 



Other surveys were carried out by Banff et al. ( 1993 ) 
with a sample of 22 spirals located mainly in the virgo clus- 
ter and by Delgado & Perez (1997) with 27 spiral galax- 
ies. In the later survey 25% of objects were identified as 
having double slope LFs. More detailed studies have also 
been performed on single objects. Scowen et al.(1992) have 
studied M101 and found an index of a = — 1.85 ± 0.05 for 



a power law fit to the LF. Rand (1992) examined M51, 
Walterbos & Braun J1992D M31 and Wyder et al. ( [1997D 
M33. Also an example of a barred galaxy, NGC 7479, was 
studied by Rozas et al. (1999), where evidence for a dif- 



ference in the properties of the regions located in the bar 
from that of the overall disk were found. 

Studies of the Ha emission in grand design spirals were 
carried out by Cepa & Beckman ( 1989| , 1990| ), Knapen 
et al. ( 1993| ) and continued in a series of four grand de- 
sign spirals by Rozas et al. ( [19964 |l996b| , [1998^ . The 
LFs in all four cases showed a clear slope transition at 
logifja = 38.6erg/s accompanied by a local maximum 
or "glitch" . In addition to the LF break, they also found 
a break at the same luminosity in the surface brightness 
gradient of the Hn regions. The gradients of less lumi- 
nous regions were about constant, whereas the gradients 
of more luminous regions were steeper and growing with 
luminosity. Yet another break has been reported at this lu- 
minosity in the central surface brightness (Beckman, priv. 
comm.). T hese f eatures have been interpreted by Beck- 
man et al. ( 1998 ) as representing a physical change in the 
Hn regions at this characteristic luminosity. Specifically, 
they contend that a transition from ionization to density 
bounding occurs at logina = 38.6erg/s. They have also 
proposed the use of this feature as a standard candle. 

Related studies of the diffuse ionized interstellar 
medium (DIG) have been carried out by Walterbos & 
Braun (1994), Ferguson et al. ( |1996[ ) and Greenawalt et 
al. (1997). The diffuse emission was found to be spatially 
associated with the Hn regions and exhibited a variety 
of different morphologies. The fraction of emission com- 
ing from this diffuse component was estimated at 20- 
50% (Ferguson et al. 1996). Leaking photons from Hil 



This paper presents statistics on the Hn regions in 
M101 from HST observations. The very high resolution of 
the images will allow us to investigate in detail the role 
played by resolution. We will find that the linear resolu- 
tion required to perform reliable measurements of central 
surface brightness, radial gradients, and diameters is quite 
stringent and can be accomplished from the ground with- 
out the aid of adaptive optics for only the nearest galaxies. 
The effects of blending on the luminosity function will also 
be found to be significant, with the diffuse emission mak- 
ing an important contribution. 

Also examined is the intrinsic clustering characteris- 
tics of the regions and diffuse emission and the role this 
plays in shaping the behavior of the blending. Evidence for 
a possible break at a characteristic clustering regime will 
be presented, which when combined with blending effects 
may provide a mechanism for producing the LF disconti- 
nuity along with it's relative stability in luminosity. 

In Sec. | of this paper we present the data and reduc- 
tion methods. A description of the methods used to mea- 
sure the properties of the H n regions is given in Sec. ||. The 
results and comparisons to some simple analytic models 
are then given in Sec. [|, followed in Sec. || by a discussion 
of their implications. 

2. Observations 

The galaxy studied here is the nearby grand design Sc 
galaxy M101. It was chosen because it is nearly face on 
and lies at the relatively close distance of about 7.4 Mpc 



(Kelson et al. 1996). These properties limit problems due 
to projection effects and provide for very high linear spa- 
tial resolution. 

The observations of M101 were obtained with the 
HST-WFPC2 on the 15th and 16th of September 1994 
and on April 4th 1996. Both frames are public and have 
been retrieved from the HST-archive. Fig. [I] shows their 
approximate location in the outer spiral arms. The ground 
based image was taken by Sandage ( 196l| ) in blue light 
(3300-5000 A) at Mt. Palomar. 

The pixel size is O'.'l for the WF-Chips and defines 
the effective resolution, i.e. the PSF FWHM is smaller. 
The Ha (F656N non-shifted) data consists of two sets of 
two exposures. The sets have combined exposure times of 
4000 s and 3600 s. Continuum subtraction was carried out 
either with exposures taken with the wide R-Filter F675W 
(2x2000 s) or with the medium bandwidth V-Filter F547M 
(2 x 400 s). 

Besides the HST pipeline calibration, we used the 
data reduction package IRAF and especially the package 
STSDAS for further processing. Particular care was re- 
quired to remove the cosmic rays through a multiple step 
process. All frames were normalized to a count flux of 
10~ 18 ergs -1 cm" 2 using the method described in Holtz- 



regions, turbulent mixing and shock ionization have been 
suggested as possible ionization mechanisms. 



man et al. ( 1995 ). Absolute fluxes were calculated using a 
distance of 7.4 Mpc. 



Fig. 1. Ground based frame of M101 with the location of the two WFPC frames. 



The linear scale of the HST observations with a pixel 
scale of O'.'l and distance of 7.4 Mpc is 3.6pc/pixel. For 
comparison with typical ground based imaging conditions, 
we also generated a set of re-pixeled frames based on a 
distance of 20 Mpc and pixel scale of 0'.'28. This gives a 
linear scale of 27pc/pixel and corresponds to a magnifi- 
cation factor of 7.5. The frames were further processed by 
convolving with a gaussian corresponding to a seeing of 
0'.'8 FWHM, leading to a resulting linear scale of 77.6 pc 
FWHM. 

In Fig. ^| is an example of a high resolution Ha and 
V-Band image, along with several additional Ha frames 
which have been shifted to varying distances using the 
above ground based seeing conditions of 0'.'8 FWHM. The 
top left and right panels correspond to the high resolution 
Ha and V frames, respectively. The lower four panels cor- 
respond, in clockwise order, to distances of 10, 15, 20, and 
35 Mpc and linear scales of 38.8 pc, 58.2 pc, 77.6 pc and 
135.8 pc FWHM, respectively. 



3. Methods 



3.1. Till-Region Definition, High Resolution 

As the first selection criterion for Hn regions we adopt 
the requirement that they consist of at least nine contigu- 
ous pixels, each with an intensity of at least three times 
the r.m.s. noise level over the local background, which in- 
cludes any contribution from local diffuse emission. With 
an overall r.m.s. of about 2 normalized counts, the mini- 
mum luminosity according to this definition is 10 35 6 erg/s 
and the minimum diameter is 10.8 pc. 

At high resolution many blobs of diffuse Ha emission 
with various shapes can be seen which at lower resolution 
are no longer clearly discernable. These "false" regions 
must be discriminated from the "true" Hil regions. We 
therefore adopt a second criterion that every low surface 
brightness region with clearly unsymmetric morphology 
must have a central stellar counterpart in the V band. If 
in a multi-peaked complex we can detect more than one 
well separated ionizing star cluster, it is separated into 
multiple regions, otherwise the complex is treated as a 
single entity. 

Small round high surface brightness regions without 
detectable ionizing star cluster are included in the HII re- 
gion definition, because we could not exclude the possibil- 
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Fig. 2. A visual illustration of the effects of resolution on the Hn region population. The two upper panels show 
the same area in a high resolution HST Ha (a) and V band (b) image with 3.6 pc pixelsize. The bottom four panels 
correspond, in clockwise order, to an equivalent ground based Ha image at distances of c 10, d 15, e 20, and f 35Mpc 
with an adopted seeing of ft' 8 FWHM. The resulting linear scales are c 38.8 pc, d 58.2 pc, e 77.6 pc and f 135.8 pc 
FWHM. 
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ity that a central cluster lies below our V band detection 
limit or is hidden by dust. This definition is more precise 
than the one of Rand Q1992Q and Rozas et al. ( |1996a| ), 
who define every peak of emission as an Hii region. Re- 
gions that meet the first criterion, but fail the second one 
are marked as diffuse regions. Still, these account for only 
8% of the detected diffuse emission. 



3.2. Hii Region Definition, Low Resolution 

For the smeared frames we adopt the same first condi- 
tion as in the high resolution case regarding size and flux. 
Therefore, with an overall r.m.s of about 15 normalized 
counts the minimum luminosity is 10 36 ' 4 erg/s and the 
minimum diameter 81 pc. This time we count every emis- 
sion peak as an H II region. 

3.3. The REGION Software 

To find the statistical properties of the H n re gions we use 

Heller 



Data 



Do [pc 



No 



1999 



the REGIO N software package (Rozas et al. 
et al. |l999|) . This semi- automated technique identifies the 
H II regions, determines their position, and measures their 
area and flux, corrected for the local background. Further 
processing can be applied manually: deleting, adding, sep- 
arating, and the combining and rounding of regions that 
do not meet our definition. 

We find 237 and 162 Hn regions, respectively, in the 
high and low (smeared) resolution frames. 

4. Results 

4-.1. Overall Smearing Effects 

A first impression of the effect of resolution can be ob- 
tained from Fig. |^, where differences between the high 
resolution frames and the smeared ones in terms of the 
H II region sizes and morphologies are clearly visible. 

A first quantitative approach is taken by tabulating the 
effects that the smearing has on the identified population: 

1. An Hii region and diffuse emission are smeared to- 
gether to form a more luminous region (31%). 

2. More than one Hn region and diffuse emission are 
smeared together to form a more luminous region 
(28%). 

3. Diffuse emission is reshaped by the smearing and de- 
tected as an Hii region (41%). 

4. Small Hii regions disappear in the noise. 

The percentages refer to the number of identified regions 
in the smeared images. It should be noted that more than 
2/3 of the regions in the smeared case actually consist of 
either multiple regions blended together or are not genuine 
H II regions, in that they lack an internal ionizing source. 

We also observed in the high resolution frames a class 
of compact high surface brightness regions which show 



a low resolution 64.8 1285 

b low resolution corrected 70.0 687 

c high resolution 29.2 318 

Table 1. Parameters of the integral diameter distribution 



a very high degree of symmetry. They are often found 
close to a large irregular region or complex. Because of 
their distinct morphologies these were considered true H n 
regions in the count, even if an internal ionizing source was 
not observed. 

4-. 2. Diameter Distribution 

The integral diameter distribution shown in Fig. ^ de- 
monstrates the clear difference between the resolved and 
smeared frames. The top panel shows the smeared case, 
the middle panel the smeared case excluding those re- 
gions which are dominated (> 80%) by diffuse emission, 
and the bottom panel the resolved case. Whereas the re- 
solved regions show a size distribution from about 10 pc 
to 220 pc, the range in the smeared case includes regions 
from roughly 150 pc to greater than 400 pc. Parameters 
determined by fitting to the exponential form, 



N = N c D « . 



(1) 



(Van den Bergh 1981) are given in Table 1. 

A plateau-like structure is seen in the resolved distri- 
bution at the high luminosity end. In the high resolution 



study of Wyder et al. ( 1997) a similar distribution is found 



for M33, as well as a similar range of diameters. However, 
no obvious corresponding feature is seen in the plot of 
log£) vs. log!/ given in Fig. [|. The resolved and smeared 
data can be fit by a power law with indices of 0.39 and 
0.25, respectively. Substituting this relationship into the 



power law luminosity function of Sec. 4.3 gives a power law 
diameter distribution with an expected index of f3 — —2.9. 
An earlier usage of this functional form for the di amete r 
distribution can be found in Kennicutt & Hodge ( 1980 ). 
Fig. |i| demonstrates that the two distributions only ap- 
proach each other at rather high luminosities, indicating 
that the measured diameters of regions of all but the high- 
est luminosities are significantly modified at ground based 
resolutions. 

In Fig. |5| we show the corresponding differential diam- 
eter distribution along with a power law fit of the form, 

dN(D) = BD P dD, (2) 

with index (3 — — 2.84±0.16. The resolved data are clearly 
better represented by the power law than by the exponen- 
tial form. This is in contrast to the lower resolution data, 
where the exponential is a better representation. There- 
fore resolution effects do not only dramatically change the 
characteristic diameter D$, but can affect the apparent 
functional form of the distribution as well. 
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Fig. 3. Integral diameter distribution. Panel a shows 
the diameter distribution for the smeared frames (filled 
squares, linear resolution 77.6 pc FWHM) compared with 
the resolved distribution (small open circles, linear res- 
olution 3.6 pc pixelsize). Panel b shows the distribution 
for the smeared frames excluding regions dominated by 
diffuse emission. Panel c gives the distribution for the re- 
solved frames. Exponential fits are shown with a solid line. 



4-3. Luminosity Function 

The combined luminosity function (LF) for the full set of 
resolved frames is given in Fig. ^|c. It displays the char- 
acteristic power law form along with a flattening at lower 
luminosities. The flattening may be explained as due to a 
transition from ionizin g stel lar clusters to single ionizing 
stars (Kennicutt et al. 1989 ). It should be noted that due 
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Fig. 4. Logarithmic plot of region diameters against lu- 
minosity in the smeared (solid circles, linear resolution 
77.6 pc FWHM) and resolved (open circles, linear resolu- 
tion 3.6 pc pixelsize) cases. Both distributions can be fit by 
a straight line which represents a power law. The measured 
diameters of regions of all but the highest luminosities are 
significantly modified at ground based resolutions. 
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Fig. 5. Differential diameter distribution along with a 
power law fit for the high resolution frames with 3.6 pc 
pixelsize 



to completness effects most ground based studies find a 
turn-over instead of a flattening. 

If we include the diffuse regions in the LF, we obtain 
a turnover at lower luminosities rather than a flattening 
although the power law index changes little. In particular, 
since there are no diffuse regions with logL > 38.0 erg/s, 
the high luminosity part remains unaffected. 

Parameters of a power law fit (Kennicutt et al. 1989| ), 



diV(L) = AL a dL, 



(3) 
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Fig. 6. Differential luminosity functions. Panel a shows 
the LF of the smeared frames (full squares) in compari- 
son with the resolved function (open circles) . Panel b dis- 
plays the LF of the smeared frames excluding the diffuse 
emission dominated regions. Panel c gives the LF of the 
resolved frames along with a power law fit. 



are given in Table 2. The LF is consistent with a single 
slope with exponent a = —1.74 ± 0.08 in the luminosity 
range logi = 36.7 to 39.3 erg/s. 

The corresponding LF for the smeared frames are given 
in Fig. |6|a and b. Unlike the resolved frames, the LF for the 
smeared frames may be fit consistently with either a single 
or double slope. Unfortunately the lack of better statistics 
does not allow us to make a clear determination of whether 
a break in the LF does actually occur in this case. However 
two effects can be clearly distinguished: a decrease in the 



a low resolution -1.7 37.2 

b low resolution corrected -1.6 37.3 
c high resolution -1.74 36.8 

Table 2. Parameters of a power law fit to the luminosity 
functions 
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Fig. 7. Logarithm of the ratio of the number of regions in 
the smeared to resolved case for a given luminosity bin. 



number of regions in the lower luminosity bins and an 
increase of regions in the intermediate luminosity range. 
This trend becomes even clearer when regions dominated 
by diffuse emission (>80% in L) are excluded. 

The presence of these two resolution effects is more 
clearly demonstrated in Fig. [?], where we plot the ratio of 
the number of regions in the smeared to resolved case for a 
given luminosity bin. In this way one can easily see where 
there is an excess or deficit of regions as compared to the 
resolved numbers. The figure shows a deficit at logL < 
37.6 and an excess between \ogL — 37.6 and 38.6, which 
peaks at logi w 38.4. At higher luminosities the resolved 
and smeared LFs appear to converge, at least within the 
statistical uncertainties. 

These differences, combined with the additional effects 



discussed in Sec. |4.1| , provide a picture of what happens 
as the resolution is degraded. In general, regions are com- 
bined with diffuse emission and other regions, shifting 
them to higher luminosities. While this process has been 
discussed by other authors, the degree to which this oc- 
curs at typical ground based resolutions has not been fully 
appreciated. 

Furthermore, we propose that a break or significant 
reduction in this process occurs at luminosities above 
\ogL 38.6 erg/s. This point is demonstrated in Fig. |^ 
where we have plotted the luminosity of regions identi- 
fied in the smeared case against their brightest member 
as measured from the high resolution frames. Many re- 
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resolved log L [erg/s] 

Fig. 8. Region luminosity in the smeared case vs. bright- 
est member from resolved case. Symbols represent diffuse 
dominated regions (triangle) , single regions combined with 
diffuse emission (dots) , and blended regions together with 
diffuse emission (star). 



gions are shifted towards a limit of logi = 38.5 erg/s due 
to blending with other regions and diffuse emission. This 
mechanism may account for a break in the LF often ob- 
served around this limiting luminosity. 

4-4- Radial Profiles 

At high resolution the radial profiles typically display 
a core-wing structure. The profiles in the smeared case 
appear gaussian, with some core-wing structure in the 
brightest regions. 

4-5. Surface Brightness Gradient 

To compare the surface brightness gradients of the re- 
solved regions with earlier observations (Rozas et al. 



199S), we need to determine a single gradient based on the 



profile. We decided to fit a single gaussian to the core/wing 
profile and define the gradient as the slope at the half max- 
imum radius. 

We adopt the same definition for the smeared regions, 
where it is even more appropriate. We also tried fitting 
a linear function and found that the difference compared 
with the gaussian is about 20%, which is much less than 
the intrinsic dispersion. 

The surface brightness gradients of the resolved regions 
are shown in Fig. ^|. In contrast to the ground based work 
this distribution does not show a bilinear structure with a 
break at logi = 38.5 erg/s. Rather it exhibits a constant 
upper limit at all luminosities with a lot of scatter toward 
steeper gradients. No other clear relation is visible. 
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Fig. 9. Surface brightness gradients of the resolved Hn 
regions. 



However, the number of regions with luminosities 
greater than where the break would be is small, making 
the identification of a possible break difficult. Because of 
this we bolster our position with a simple analytic model 
of the effect of resolution on the gradients. The model 
consists of 2-dimensional convolution between two gaus- 
sians, one representing a resolved profile Fq, the other the 
point-spread function (PSF) of the effective resolution or 
seeing, 
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A reasonable assumption for the model is to take the 
value for the gradient in the resolved case, Go? &s a con- 
stant denned by the upper band of points in Fig. § Solving 
for of in Eq. g gives, 



3 v / 21n2 T 

<Jr = : — — L :- 



4tt Go 
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(7) 



Substituting this into the smeared profile, Eq. ^, we obtain 
the luminosity dependent expression, 
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Fig. 10. Surface brightness gradients. The dashed line 
represents the best fit for the constant gradient of the 
resolved regions. The solid line is the theoretical curve for 
the smeared regions. The symbols represent measured gra- 
dients from the resolved (circles) and smeared (triangles) 
frames. 
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Fig. 11. Theoretical curves showing the resolution effects 
on the measured surface brightness gradients. The dashed 
line represents our constant best fit to the resolved data. 
Solid lines show the expected gradients under ground 
based conditions for distances between 5 and 35 Mpc, cor- 
responding to linear resolutions between 19.4 pc FWHM 
and 135.8 pc FWHM. The steps in distance are 5 Mpc, 
increasing upwards. 



We adopt a model resolution based on good ground 
seeing conditions of FWHM = 0'.'8 and a distance of D = 
20 Mpc, giving a value of a s = 32.9 pc. To find the constant 
fc r , we fit the relation between measured of and luminos- 
ity, obtaining a value of k T = 3.09 x 10 - 34 pc 3 erg _1 s. 

The results of these calculations are plotted in Fig. |l^, 
together with the measured gradients for both resolutions. 
While the statistics are lacking, the overall trend of the 
gradients in the smeared case follows that found by Rozas 
et al. fll99l| ). They argue that this is evidence for a transi- 
tion in the physical properties of the Hil regions. However, 
the simple model shows that the gradients measured in the 
smeared frames are in fact completely dominated by the 
effective seeing at the lower luminosities and do not rep- 
resent the intrinsic properties of the clouds. The break in 
the gradients which is seen in lower resolution observations 
marks the zone of transition where the size of the regions 
is of the same order as the effective seeing gaussian. The 
fact that the measured gradients lie above the theoretical 
curve is likely due to the neglect of the contribution from 
diffuse emission and blending in our simple model. 

In Fig. ^ we plot various model curves representing 
the resolution effects on the measured surface brightness 
gradients for different distances at a seeing FWHM of 
0'.'8. For the given range of luminosities, all galaxies at 
distances larger than about 10 Mpc display seeing domi- 
nated radial profiles, corresponding to a linear resolution 
(FWHM) of around 40 pc. 



4--6. Central Surface Brightness 

In this section we consider the central surface brightness 
S c of the Hn regions. Beckman et al. (priv. comm.) find 
that for regions below logi = 38.6, S c is proportional to 
L 1 / 3 , while for regions above this value the central surface 
brightness rises more steeply and shows a larger scatter. 

For a theoretical treatment similar to what we have 
done with the gradients, a more reliable estimate of the 
core profile is required. If we use our previous fits to the 
entire region we would greatly underestimate the central 
brightness. Therefore we determine S c from a gaussian 
profile that is fit only to the core area of each region in 
both the resolved and smeared frames. 

The core profile may be defined as 



F CI (r) 



L 



2-KOr 



(9) 



Next we convolve this core profile with the seeing gaussian 
as defined in Sec. 4.5, 



F cs (r) 



L 



27T (<7 c 2 + CT S 2 ) 

Substituting the approximate relationship 



(10) 



— k r L 



(11) 
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Fig. 12. Central surface brightness. The dashed line rep- 
resents the best fit of the resolved data. The solid line 
shows the theoretical curve for the smeared regions. Open 
circles are the resolved data points, solid circles represent 
the smeared data. 



into Eqs. || and [T(], we find 
L L3 



2na r 



27TfcJ 



and 



L 



27r(a c 2 +a s 2 ) 2lT (fc ci )l +(Ts 2 



(12) 



(13) 



As determined from the resolved data, we take a value for 
k c of e.SxlCT^p^erg^s. 

The results of the calculations are plotted against L 1 / 3 
in Fig. [l^, together with the data for both resolutions. 

The two datasets exhibit large differences. The 
smeared data show a nearly constant behavior with lit- 
tle scatter up to log!/ rj 38.6, beyond which S c rises 
more steeply and shows larger scatter. This behavior 
closely matches that of ground based data (Beckman, priv. 
comm.). In contrast, the resolved data exhibit a larger 
scatter at all luminosities. 

Considering the rough approximations used in deriving 
the theoretical relationship for the smeared regions, the fit 
to the smeared data is remarkably good. The correspon- 
dence could be further improved if the effects of blending 
and diffuse emission were taken into account. 

This simple model provides a reasonable description 
of the characteristics of the central surface brightness in 
the lower resolution ground based observations. It demon- 
strates that resolution effects are important not only for 
measurements of gradients, but also for the determination 
of central surface brightness. 

In Fig. |l3j we plot various model curves represent- 
ing the resolution effects on the measured central surface 
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Fig. 13. Theoretical curves showing the resolution ef- 
fects on the measured central surface brightness. The 
dashed line represents a possible fit to the resolved data. 
Solid lines show the expected central surface brightness 
under ground based conditions for distances between 5 
and 35 Mpc, corresponding to linear resolutions between 
19.4 pc FWHM and 135.8 pc FWHM. The steps in dis- 
tance are 5 pc, decreasing upwards. The transition from 
a nearly constant value to increasing brightness is clearly 
visible. 



brightness for different distances at a seeing FWHM of 
0'.'8. Again we find that for galaxies at distances greater 
than about 10 Mpc the central surface brightness is domi- 
nated by the seeing. This corresponds also to a linear scale 
of about 40 pc (FWHM). 

4-7. Minimal Spanning Tree 

In order to investigate the clustering properties of the re- 
gions we compute a Minimal Spanning Tree (MST) based 
on the spatial locations of the Hll and diffuse regions 
(Sec. |3.2p . The MST for any given set of points consists of 
a unique set of edges connecting the points, such that the 
sum of the edge-lengths is a minimum. The MST provides 
a reproducible and unbiased graph that contains informa- 
tion on the intrinsic linear associations in the point set. 
For the analysis of linear features the MST has been found 
to be superior to other methods, such as the two point cor- 



relation function (Barrow et al. 1985). 

We have chosen to include the diffuse regions in the 
analysis because they are a good tracer of the spatial dis- 
tribution of the diffuse emission in general, even though 
they represent only 8% of the total diffuse flux. A graph- 
ical depiction of the region spatial distribution and the 
corresponding MST are shown in Fig. |TJ for the two dif- 
ferent WFPC2 fields in the galaxy. 

The spatial distributions exhibit some clear features. 
The regions are not located uniformly over the frames, 
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Fig. 14. a and b Region spatial distribution for the two fields. The different symbols depict the H n regions (solid 
circles), diffuse regions (open circles), and high luminosity Hn regions with Ljja > 10 38 erg/s (dotted circles), c and 
d Corresponding MSTs. Empty region at bottom right of each panel represents location of PC chip. 



but rather lie along ring- or bubble-like structures. Fur- 
thermore, on the first frame four particularly luminous 
regions lie along a common edge of two rings. An expla- 
nation for this distribution could be induced star forma- 
tion in an expanding supernovae shell. A rough estimate 
of the radii of these structures is 0.7-1.0 kpc. This is of 
the order of the largest HI holes observed in gala xies, the 



so-call ed "su perbubbles" (Kamphuis et al. |1991[ Oey & 
Clarke 1997 ). Another possible formation mechanism for 
the observed features might be through gravitational and 



thermal instabilities in a differentially rotating disk (Wada 
& Norman |1999[). 



Similar distribution features are also observed for HI 
gas. A study of 11 nearby spirals by Braun (1997) has 



found that between 60-90% of the total HI line flux comes 
from a network of filamentary features with typical widths 
of about 150 pc, which is comparable to the region distri- 
bution in Fig. [lj. 
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Fig. 15. Integral MST edge-length distribution of the two 
fields combined at a linear pixel scale of 3.6 pc. 



4.7.1. Resolution effects 

To model the resolution effects on H n regions, we adopt a 



Gaussian profile for the regions similar to Sec. 4.5. For sim- 



plicity of the calculations we discuss the effect on two iden- 
tical region profiles. A simple calculation shows that the 
peaks of two identical gaussian profiles remain resolvable 
under convolution with a seeing gaussian if the peaks are 
originally separated by a distance larger than 2y 'of + of. 
Only for the most luminous regions is cr r greater than a s . 
Because the luminous regions are not clustered close to- 
gether (see Fig. |l4|), we can as a resonable approximation 
drop <7 r . Therefore regions that are closer than two see- 
ing er s can not be separated. The pixelsize places further 
constraints on the resolvability of regions. In order to dis- 
criminate two emission peaks from each other, they must 
be separated by at least one pixel of lower flux. Hence the 
separation between the peaks must be at least two pix- 
els. Our adopted ground based conditions are given by a 
pixelsize and seeing FWHM of respectively, 0'.'28 and 0'.'8, 
which correspond to a s = f/.'34. Then the limiting quan- 
tity is 2er s of 0'.'68, corresponding to a linear dimension of 
about 66 pc at a distance of 20Mpc. 

The integral edge-length distribution, shown in Fig. |l5|, 
gives directly the fraction of edges which have lengths 
above a given value. We can see that about 45-50% of 
all edges cannot be resolved at our adopted ground based 
resolution, demonstrating the importance of blending. The 
degree of blending will depend not only on the resolution, 
but also on the clustering properties of the Hil regions 
and diffuse emission. 

In Fig. |l^ we present the differential MST edge-length 
distribution along with power law fits. The parameters for 
the fits are given in Table 3. The power law index varies 
little between the frames (1.5-1.7), but the critical length 
where the distribution flattens differs substantially (16- 
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Fig. 16. Differential edge-length distribution for the MST 
of H II and diffuse regions for the individual fields (a & b) 
and for both fields combined (c) at a linear pixel scale of 
3.6pc. 



field 



Slope log l m in lmin [pc] 



a first 
b second 
c combined 



-1.49 1.20 16 

-1.69 1.65 45 

-1.53 1.40 25 

Table 3. Parameters of power law fit to the edge-length 
distributions. 



45 pc). This could be a hint for a different minimum clus- 
tering regime for the two regions in the galaxy. This range 
of values is below the resolution of 66 pc of our smeared 
frames. 
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Fig. 17. Number of clusters vs. clipping distance. The 
symbols represent the clusters from the individual fields 
(open circles) and the fields combined (solid circles). Also 
given along the top is the equivalent galaxy distance. 

To further quantify the impact of region clustering 
combined with insufficient resolution, we simulate the 
blending of regions by removing the edges that are longer 
than twice the seeing er s . The clipping length l c is now 
treated as a function of the distance D, 



lc 



7T • 0"68 
3600 • 1 80 



D. 



(14) 



We define a cluster as two or more regions which remain 
connected by edges following the clipping. With this defi- 
nition all the regions in a given cluster would be smeared 
together into a single region at the specified seeing and 
distance. 

Fig. |l7| shows the number of clusters as a function of 
the clipping distance and the equivalent galaxy distance, 
with an assumed seeing of 0'.'8. The curves first increase 
rapidly, reach a maximum number of clusters at around a 
clipping length of 40 pc, and then slowly decrease as the 
clipping becomes larger. In Fig. [l8]the corresponding aver- 
age number of regions per cluster is given. While differing 
in detail, the overall clustering characteristics of the two 
fields are similar, growing approximately linearly with the 
clipping (smearing). The difference in slopes could be sta- 
tistical, but also may have a physical origin in differing 
clustering properties of the two fields. 

The clustering behavior may be interpreted in terms of 
two physical spatial regimes. Regions that are separated 
by edges lying in the first regime of 40 pc are thrown to- 
gether by the smearing, creating a large number of clus- 
ters. This first regime corresponds to the flat part of the 
differential edge-length distribution. With increasing cut- 



clipping [pc] 

Fig. 18. Average number of regions per cluster vs. clip- 
ping distance. The symbols represent the clusters from 
the individual fields (open circles) and the fields combined 
(solid circles). Also given along the top is the equivalent 
galaxy distance. 

off length the clusters, which are now grouped in the sec- 
ond regime, are slowly blended together according to their 
power law distribution. This behavior can been seen by eye 
in Fig. |. 

A check may be made by determining the sum of 
the numbers of clusters and unclustered regions with 
logi > 36.7erg/s at a clipping of 60pc. This sum should 
be roughly equal to the number of regions in the smeared 
frames. This is because a clipping of 60 pc corresponds 
to a distance of approximately 20 Mpc (Eqn. ^J) and the 
smallest detected regions in the smeared case have lumi- 
nosities of logL = 36.7erg/s. The totals are in fact quite 
close, with values of 162 and 157, respectively. 

If the blending of regions were responsible for the ob- 
served break in the luminosity function, with these cluster- 
ing properties one would expect that the break would first 
appear at distances of around 10 Mpc (Of! 8 FWHM) and 
then only change slowly with increasing distances. The 
corresponding linear scale again is about 40 pc FWHM. 

4.7.2. Physical Clustering 

Putting resolution effects aside, we now examine the phys- 
ical spatial clustering properties of the Hn and diffuse 
regions. In Fig. |l9| we show the separate differential edge- 
length distributions for the H n regions and diffuse regions. 
The slopes of the power law fits are a = — 1.17 ± 0.07 and 
a = —1.29 ± 0.12, respectively. The flattening of the dis- 
tributions occurs at edge-lengths of about 20 and 32 pc. 
These values do not differ substantially, indicating a sim- 
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5. Discussion 
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Fig. 19. Differential MST edge-length distribution of a 
H II and b diffuse regions from both fields combined at a 
linear pixel scale of 3.6 pc. 



ilarity in the clustering properties of the Hn and diffuse 
regions. 



4.8. The Ha Flux 

To make a quantitative estimate of the diffuse emission 
contribution, we subtract the contribution of the Hn re- 
gions from the total Ha flux. Since the local background 
has been removed from the region flux, our measurement 
of the region contribution corresponds to a lower limit. To 
obtain an upper limit we subtract the region background 
fluxes as well. We find that 35-39% of the Ha flux in the 
resolved frames is attributable to individual H II regions. 
For comparison, Ferguson et al. ( |1996[ ) find values of about 
50-70% and Rozas et al. ( |1999[ ) 40-65%. The remaining 
flux must come from the diffuse regions and background. 
If we determine the same ratio for the smeared frames we 
obtain a somewhat different result. Here 64-74% of the 
total emission is found to be attributable to Hn regions 
and only 26-36% to the diffuse part. 

A substantial fraction of the total diffuse emission is 
included in the Hn region luminosities by the effect of 
smearing. This is another indication of the importance of 
resolution in the interpretation of H II region data. 



We have shown in the previous sections that the basic 
properties of H II regions determined from high resolu- 
tion HST frames differ significantly from even high quality 
ground based measurements. Only a few previous studies 
have had the necessary linear resolution to make reliable 
measurements of their statistical properties, such as size 



and luminosi ty di stributions (Walterbos & Braun |1992 



Wyder et al. |l997p 



We simulate the effects of lower resolution on our HST 
frames by artificially smearing the data to typical ground 
based resolutions. The diameter distribution from high 
resolution HST observations exhibits a totally different 
range of values compared to those determined from lower 
resolution data. We find a diameter distribution which is 
more consistent with a power law than with the often used 
exponential function. This puts into doubt the concept of 
a characteristic size for Hn regions. 

The luminosity function can be fit by a single slope 
power law above a luminosity of logi = 36.8erg/s. Below 
this value the distribution is flat. Unfortunately, the exis- 
tence of a break at higher luminosities can not be reliably 
determined through direct examination of the luminosity 
function derived from the M101 data because of the small 
number of regions with luminosities above the expected 
location of the break. In order to make such a determina- 
tion, greater coverage of M101 will be needed to increase 
the statistics, particularly at the high luminosity end. 

The LF of the artificially smeared data (Fig. ||a) is 
flatter in the mid-luminosity range. A clear identification 
of a high luminosity break can also not be reliably estab- 
lished. A close comparison of the frames at the two dif- 
ferent resolutions has given us important insight into the 
influence of resolution on the measured H n region prop- 
erties. We find that as the resolution is reduced regions 
which are spatially close together are blended into a single 
observed "region" , along with a significant amount of dif- 
fuse emission. This process leads to the shifting of regions 
from low luminosity bins to higher bins. A luminosity of 
logL w 38.5erg/s appears to mark an upper limit for this 
luminosity shifting effect. The blending of regions alone is 
not sufficient to explain this limiting luminosity, but re- 
quires the inclusion of the diffuse emission. While still not 
fully understood, it appears that the intrinsic properties 
of the diffuse emission clouds, such as their size, surface 
brightness, and spatial association with the Hn regions, 
play an important role. Modeling of the LF with the inclu- 
sion of the various effects of degraded resolution (Sec. 4.1 ) 
would be helpful in determining the nature of the high 
luminosity break and its possible relation to the diffuse 
emission. 

The idea that the high luminosity break is a result 
of the blending of Hn regions and the masquerading of 
diffuse regions as genuine ionization sources differs from 
the hypothsis that it is due to a transition in the physical 
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state of the Hn regions from ionization bounded to den- 
sity bounded (Rozas et al. 1996a; Beckman et al. 1998j ) . 
Still a connection might be given by the nature of the 
diffuse emission. If the diffuse emission close to the Hn 
regions is only ionized by the nearest source, then it may 
be considered as being associated with the regions. When 
a region and its associated diffuse emission are thrown to- 
gether by the effects of resolution, then the blended region 
may be seen as being ionization bounded when sufficient 
diffuse gas is included. On the other hand a region may be 
seen as density bounded when the amount of diffuse gas 
included is insufficient with respect to the luminosity of 
the internal ionizing source. In this scenario the amount 
of diffuse gas included in a blended region is controlled by 
the resolution and hence to some extent so would be the 
position of the high luminosity break in the H n region 
LF. 

With even good ground based data the definition of 
what is an H n region is heavily influenced by the reso- 
lution. At high resolution one can see that these smeared 
"regions" really consist of a complex of associated struc- 
tures, such as diffuse clouds, arclets, shells, and filaments. 
Hence the surrounding gas is far from being evenly dis- 
tributed around the ionizing source. With such a non- 
uniform topology it would be difficult to achieve a true ion- 
ization bounding, because the possibilities for a Lyman- 
continuum photon to escape unhindered are numerous. 
Such a situation is a poor representation of an idealized 
Stromgren sphere. It is also not clear to what degree exter- 
nal sources may contribute to the ionization of the diffuse 
component and particularly of the diffuse regions. 

A direct examination of a sequence of artificially reso- 
lution degraded frames, and confirmed by the MST anal- 
ysis, shows that the effects of smearing become particu- 
larly prominent around a critical linear resolution of about 
40 pc FWHM and thereafter only worsens slowly. If the 
high luminosity break in the Hn region luminosity func- 
tion is an artifact of resolution then this characteristic be- 
havior may explain the stability of the break luminosity 
over a relatively large range of galaxy distances. Therefore 
the stability of the break luminosity is not unambiguous 
evidence against blending as a mechanism for producing 
the discontinuity. Clearly, modeling of this mechanism will 
be required to determine its viability. However, the impor- 
tance of the diffuse emission should not be ignored and any 
model of the blending will need to take this into account 
along with the characteristics of the clustering. 

The underlying reason for this behavior seems to be 
the clustering characteristics of the Hn regions and dif- 
fuse emission. The MST edge-length distribution exhibits 
a power law form with a flattening at about 25 pc. The 
length at which this transition occurs may be the pri- 
mary factor influencing the location of the high luminos- 
ity break. The physical reason for the flattening at this 
value of edge-length is not clear, but may be related to 
a typical maximal fragmentation length within a single 



precursor gas-cloud. The power law slope would then rep- 
resent the distribution of distances between the precursor 
gas-clouds. Another possibility could be related to two 
different mechanisms of star formation as described by 
Scowen et al. ( 1996j ). Whereas the power law edge- length 
distribution might correspond to the fragmentation of the 
ISM induced by the passage of a spiral arm density wave, 
the flat part may represent typical spatial sizes of local 
subsequent star formation. In this case the patches of dif- 
fuse emission might be interpreted as expanding HII re- 
gion remnants. If the break in the LF is interpreted as due 
to these clustering properties, its frequent appearance de- 
rived from ground-based data may suggest a similarity of 
two different regimes of star formation in M101 and other 
late-type spirals. 

There is an indication that the clustering character- 
istics between the two observed fields differ. While the 
edge- length distributions show similar slopes, the value at 
which the power law flattens out differ by about 30 pc. It is 
not clear if this difference has a physical origin or is merely 
statistical. Such a difference is not found when we compare 
the luminosity functions or diameter distributions for the 
different fields. The flattening of the luminosity function is 
found for both fields at logL = 36.8, which corresponds to 
a diameter of about 30 pc. If the difference is indeed phys- 
ical in nature, then one must explain how the properties 
of the fragmentation regimes can differ without affecting 
the other star formation properties. 

Our simple analytic model of the resolution effects on 
the profiles can successfully reproduce the general behav- 
ior of the surface brightness gradients seen at lower reso- 
lutions. While the profiles of the lower luminosity regions 
(logL ^ 38.5) are dominated by the point spread func- 
tion, the profiles of higher luminosity regions which have 
spatial sizes greater than the seeing scale are expected to 
represent the true properties of the regions. We therefore 
interpret the break in the surface brightness gradient dis- 
tribution as a resolution effect. This simple model provides 
a prediction on how the luminosity of the break changes 
with distance. 

The measured central surface brightness of the re- 
solved regions shows a large scatter and does not exhibit a 
well defined trend with luminosity. This is not true for the 
smeared regions which exhibit a nearly constant value up 
to logL = 38.5erg/s. After this the central surface bright- 
ness increases steeply with increasing luminosity. Our sim- 
ple model again reproduces this behavior fairly well and 
predicts the transition luminosity will vary with distance. 

In both cases the profile and the central brightness be- 
come seeing dominated at linear resolutions worse than 
40 pc FWHM. Therefore this linear resolution marks the 
minimum resolution required for reliable determinations of 
region characteristics. This kind of resolution is currently 
only achievable from the ground for the nearest galaxies. 
Taking this into account the observed breaks at lower res- 
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olution do not seem to be reliable evidence for a transition 
from ionization bounded to density bounded Hn regions. 

Concerning the question as to whether the high lumi- 
nosity break in the LF can be used as a standard candle as 
proposed by Beckman et al. (priv. comm.), more observa- 
tions will be required to validate the stability of the break 
luminosity with distance (resolution) and galaxy type. A 
clearer understanding of its physical origin would greatly 
assist in such an effort. In our interpretation, the break is 
a result of the clustering characteristics of the H n regions 
and diffuse emission. As seen in the MST analysis, the dif- 
ferences between the two fields we studied in MI01 may be 
an indication that significant differences in the clustering 
properties exist even within a single galaxy. Depending on 
the degree to which this were true, the utility of a global 
LF of a galaxy as a measure of its star formation proper- 
ties would be compromised. 

A further indication of variations in region clustering 
may be in th e stu dy of the barred galaxy NGC 7479 by 
Rozas et al. ( 1999 ), who find differences between the bar 
Hn region population and that of the overall disk. Fur- 
thermore this galaxy shows a different behavior in compar- 
ison to the non-barred galaxies with respect to the break 
luminosity in the central surface brightness distribution. 
This could be a hint for a possible variation of H II region 
clustering properties with galaxy type. 

Clearly additional high resolution work is needed to 
both confirm and extend these results. More sophisticated 
methods should be applied to the radial profiles in order 
to better characterize their properties. Even a rough de- 
termination of cluster age from broad band imaging would 
be useful for constraining the physical origin of the larger 
scale HI structures and their relationship to the star form- 
ing clouds. 

With respect to possible further high resolution ob- 
servations, obtaining better statistics is vital for making 
further progress on the luminosity and diameter distribu- 
tions. To generalize the results, these investigations should 
be carried out for different galaxy types. A comparison 
between a barred and non-barred galaxy would be an im- 
portant test of the stability of the high luminosity break 
as the dynamics of the two systems are quite different. 

6. Summary 

In this paper we present high resolution HST Ha data 
of M101. For a comparison with previous ground based 
observations we artificially shift these frames to an effec- 
tive distance of 20Mpc by rebinning, changing the pixel- 
size to 0'/28, and convolving with a seeing gaussian of 0'.'8 
FWHM. For both resolutions we derive the integral diam- 
eter distribution, luminosity function, surface brightness 
gradients and central surface brightnesses. 

For the integral diameter distribution at high reso- 
lution we find a characteristic diameter Dq = 29.2 pc. 
However, a power law fits the data better with a slope 



f3 = —2.84 ±0.16. The luminosity function in the resolved 
case can be fit by a single slope power law with a = —1.74 
in the luminosity range between 10 36 ' 7 and 10 39 3 erg/s. 
It flattens for lower luminosities. The surface brightness 
gradients exhibit a constant limit with a lot of scatter to 
steeper gradients. Because of large scatter it is hard to 
find an appropriate function of luminosity for the central 
surface brightness. 

Comparing the diameter distribution at the two reso- 
lutions, we find a totally different range of values. Whereas 
the diameters of the resolved regions lie between about 10 
and 220 pc, the diameters in the smeared case show values 
from about 150 pc to greater than 400 pc. We also find a 
different functional form at the higher resolution, namely 
a power law rather than an exponential function. 

The LF in the smeared case exhibits a reduction of re- 
gions in lower luminosity bins, while an excess in the range 
between 10 378 and 10 38 6 erg/s. We attribute this behav- 
ior to the intrinsic clustering properties of the H II regions 
and diffuse emission combined with resolution effects. As 
the resolution is reduced spatially associated regions are 
blended together including a significant amount of diffuse 
emission. This process leads to a shifting of the regions 
into higher luminosity bins and may account for the often 
observed break in the LF at logL = 38.6 erg/s. 

We have further investigated this effect by examining 
the clustering quantitatively with a minimal spanning tree 
(MST). There we find a power law function for the edge- 
length distribution which flattens at about 25 pc and has 
an index of about a — —1.5. 

The edge-length where the flattening occurs may be 
the primary influence on the location of the break. The 
two different spatial regimes are interpreted as two dif- 
ferent regimes of star formation in M101. If the break in 
the LF is interpreted as due to these clustering properties 
its frequent appearance derived from ground-based data 
may suggest a similarity of two different regimes of star 
formation in M101 and other late-type spirals. 

The surface brightness gradients as a function of lu- 
minosity for the smeared regions may be described by a 
bilinear fit, but a simple analytic model shows that this 
feature results completely from the degraded resolution. 
With the same mathematical approach the break of the 
central surface brightness of the smeared regions is pre- 
dicted correctly, indicating that neither break can be con- 
sidered as a real physical effect. 

The minimum linear resolution needed to reliably mea- 
sure internal quantities is found to be about 40 pc FWHM. 
Such resolutions are only currently obtainable from the 
ground for the nearest galaxies. 

In our view the hypothesis that the high luminosity 
break in the LF is due to a transition from ionization 
bounded to density bounded H n regions is still not firmly 
established and resolution effects need to be better taken 
into account when making such studies. The intrinsic clus- 
tering properties of the Hn regions and the associated 
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diffuse emission combined with typical resolutions may be 
able to produce the observed break. 

Further high resolution observations will be needed to 
confirm these results and obtain reliable statistics on the 
characteristics of Hn regions in galaxies. Only then will 
it be possible to determine whether the clustering charac- 
teristics of H II regions differ with galaxy type and envi- 
ronment, and any connection it may have with the local 
dynamics. 
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